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ABSTRACT 

Aims. The aim of this work is to characterize the stellar population between Earth and the Orion A molecular cloud 
where the well known star formation benchmark Orion Nebula Cluster (ONC) is embedded. 

Methods. We use the denser regions the Orion A cloud to block optical background light, effectively isolating the stellar 
population in front of it. We then use a multi-wavelength observational approach to characterize the cloud's foreground 
stellar population. 

Results. We find that there is a rich stellar population in front of the Orion A cloud, from B-stars to M-stars, with a 
distinct 1) spatial distribution, 2) luminosity function, and 3) velocity dispersion from the reddened population inside 
the Orion A cloud. The spatial distribution of this population peaks strongly around NGC 1980 (iota Ori) and is, in all 
likelihood, the extended stellar content of this poorly studied cluster. We infer an age of ~ 4 — 5 Myr for NGC 1980 and 
estimate a cluster population of the order of 2000 stars, which makes it one of the most massive clusters in the entire 
Orion complex. This newly found population overlaps significantly with what is currently assumed to be the ONC and 
the L1641N populations, and can make up for more than 10-20% of what is currently taken in the literature as the ONC 
population (30-60% if the Trapezium cluster is removed from consideration). What is currently taken in the literature 
as the ONC is then a mix of several intrinsically different populations, namely: 1) the youngest population, including 
the Trapezium cluster and ongoing star formation in the dense gas inside the nebula, 2) the foreground population, 
dominated by the NGC 1980 cluster, and 3) the poorly constrained population of foreground and background Galactic 
field stars. 

Conclusions. Our results support a scenario where the ONC and L1641N are not directly associated with NGC 1980, i.e., 
they are not the same population emerging from its parental cloud, but are instead distinct overlapping populations. 
The nearest massive star formation region and the template for massive star and cluster formation models is then 
substantially contaminated by the foreground stellar population of the massive NGC 1980 cluster, formed about 4-5 
Myr ago in a different, but perhaps related, event in the larger Orion star formation complex. This result calls for a 
revision of most of the observables in the benchmark ONC region (e.g., ages, age spread, cluster size, mass function, 
disk frequency, etc.). 

Key words. Stars: formation - Stars: massive - Stars: pre-main sequence — ISM: clouds 



1. Introduction 

The Orion Nebula is one of the most studied objects in 
the sky, with observational records dating about 400 years 



pc (Kim et al. 2008)). Some of the most basic observables 
of the star formation process, like, 1) s tar formation rates 
(Lada fc Lada 1995; L ada et al. 2010), 2) s tar formation 



(|Muench et al. 


2008) 


tance for astrop 


hysics 



It is an object of critical impor- 



histo r y QHhlenbrand||1997|), 3) age spreads ( |Jeffries et al 



2011 Reggiani et al. 2011[), 4) the initial mass function 

(111 



massive star formati on region to Earth, the Orio n Nebula 
Cluster (ONC) (e.g. |Johnsonlp65l [Wiiker1p72| , 
the benchmark region for massive star and cluste: 



to the substellar regime JHillenbrand fc Carpenter 2000 
Muench et al.|2002|[Da~Rio et al.|2012||Hsu et al.|2012| ), 5) 



which is 
cluster forma- 



the fr action, s ize distribution, and lifetime of circumstcllar 
disks (iHiTlenbrand et al ||1998[ |Lada et al 



tion studies. Recent distance estimates to the Orion Nebula 
using paralla x put this object at abo ut 400 pc fro m Earth 



2000 Muench 



et al.||200Tj |Vicente fc Alves||2005[) , 6) their interpIaTwrth" 



massive star s QO'dell et al. | 199 3 [), binarity (|Petr et al.|1998 



(389t^ pc ( 


Sandstrom et al. 


2007 


), 414±7 


et al. 


2007 


), 43Y±19 pc ( 


riirota et al. 


2007) 



and 4I9±fc> 



Kohl er et al.||2006P, rota tion (|Herbs t et al.|2002^ " magnetic 



fields dFeigelson et al.|2003l alirJ Ty young cluster dynamics 



(Hillcnb rand & Hartmann 



1998 Furesz et al. 2008 Tobin 



* Based on observations collected at the German-Spanish 
Astronomical Center, Calar Alto, jointly operated by the Max- 
Planck-Institut fur Astronomie Heidelberg and the Instituto de 
AstrofTsica de Andalucia (CSIC). 
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the benchmark region for theoretical and numerical models 
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Fig. 2. Coverage of the optical datasets used in this study. 
The SDSS images used in this study are represented in 
green, the CFHT/Megacam images in red, and the Calar 
Alto 1.23m CCD observations are represented in light blue. 
The contour corresponding to an integrated intensity of 
13 CO of 30 K km/s is represented in white. North is up and 
east is left. The angular scale is indicated in the lower left. 
Background photograph courtesy of Rogelio Bernal Andreo 
(DeepSkyColors.com) 



of massive and clustered star formation (Palla fc Stahler 
19991 JKlessen fe Burkert|[2000{ parke et al.||2000| |Bonnell 
et al.|2001||Bate et al.|2003[|Tan et al.|2006 



2006 |Krumholz et al.||2011|7Allison fc Goodwin||201lj 



Huff & Stahler 



It is quite remarkable that within only 1.5° of the ONC 
there are several contiguous ,and likely overlapping, groups 
of young stars (see Figure [l]) , although few studies have 
tackled the entire region as a whole. Unfortunately the 
three-dimensional arrangement of star formation regions, in 
particular massive ones, is far from simple and is essentially 
unknown given the current distance accuracy to even the 
nearest star formation regions. It is clear however that the 
Orion Nebula cluster is partially embedded in its parental 
Orion A molecular cloud which in turn is inside the large 
~ 200 pc Orion star formation complex, where groups of 
young stars with ages from a few to abou t 10 Myr are seen 



( Brown et al. 1994[ Briceno et al. 2007 ) . It has long been 



suspected that a more evolved group (subgroup Ori OB lc 
including NGC 1981 and NGC 1980, |BlaauwTT964] [Warrenl 
fc Hesser||1978 > is in the foreground of the molecular cloud 
where the younger ONC popul ation (subgroup Ori O B Id) 
is still partially embedded (see|Gomez & Lada||1998l for a 



ep 
sibl 



There are two different views on the stellar popula- 
tion inside the Orion Nebula. The first suggests that the 
core of the ONC, the Trapezium cluster, is a distinct en- 
tity from the rest of the stellar population in the nebula, 
while the second, more prevalent, suggests the Trapezium 
is instead the core region of a larger clus t er em erging from 
the Orion Nebula. Herbig fc Terndrupj ( 1986 ) performed 
one of the first CCD observations of an area centered on 
the Trapezium cluster (covering ~ 9.2' 2 ), and from the ex- 
ceptional high stellar density found they argued that the 
Trapezium cluster was a distinct entity from the surround- 
ing stellar population, including the stellar population in- 
side the Orion Nebula. An oppo site view was proposed by 



large scale analysis of the possible interplay between these 
two subgroups). 



Hillenbrand & Hartmann ([1998 ) who compared optical and 
near-infrared surveys of the ONC with virial equilibrium 
cluster models to argue that the entire ONC is likely a sin- 
gle young stellar population. 

Confirming which view is correct is critical because they 
imply different formation scenarios for the ONC, and as- 
suming the ONC is typical, different scenarios for the for- 
mation of stellar clusters in general. While the first view 
implies the bursty formation of the bulk of the stars in a 
relatively small volume of the cloud, the second, by assum- 
ing a more extended cluster, calls necessarily for a longer 
and more continuos process, allowing for measurable age 
spreads in the young population, and for substancial frac- 
tions of young stellar objects (YSOs) at all evolutionary 
phases, from Class to Class III. Observationally, the first 
view argues that the Orion OB lc subgroup is a distinct 
star formation event from the Id subgroup while the second 
and more prevalent view argues that the two subgroups are 
the same population, i.e., the Ori OB lc subgroup is sim- 
ply the more evolved stellar population emerging from the 
cloud where group Id still resides. 

If the first view prevails, i.e., if the Trapezium cluster 
and ongoing star formation in the dense gas in its surround- 
ings represent a distinct population from the rest of the 
stars in the larger ONC region, then what is normally taken 
in the literature as the ONC is likely to be a superposition 
of different stellar populations. If this is the case, then the 
basic star formation observables currently accepted for this 
benchmark region (e.g., ages, age spread, cluster size, mass 
function, disk frequency, etc.) could be compromised. 

In this paper we address this important issue by at- 
tempting to characterize the stellar populations between 
Earth and the Orion Nebula. Our approach consists of us- 
ing the Orion A cloud to block optical background light, 
effectively isolating the stellar population in front of it. 
We then use a multi-wavelength observational approach 
to characterize the cloud's and nebula's foreground pop- 
ulation. We find that there are two well defined, distinct, 
and unfortunately overlapping stellar populations: 1) a fore- 
ground, "zero" extinction population dominated by the 
poorly studied but massive NGC 1980 cluster, and 2) the 
reddened population associated with the Trapezium cluster 
and L1641N star forming regions, supporting the first view 
on the structure of the ONC as described above. This result 
calls for a revision of most of the star formation observables 
for this fiducial object. 

This paper is structured as follows. In Sect, 2 we de- 
scribe the observational data acquired for this project as 
well as the archival data used. In Sect. 3 we present the 
results of our approach, namely the identification of the 
two foreground populations and its characterization. We 
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Fig. 1. Optical image of the North end of the Orion A molecular cloud, including the relatively more evolved populations 
of NGC 1981, NGC 1977 and NGC 1980 (Orion OB lc subgroup) and the Orion Nebula Cluster (Orion OB Id subgroup), 
projected against the Orion Nebula (M42). This image illustrates well the complicated distribution of young stars in the 
vicinity of the ONC, with scattered groups of more evolved blue massive stars projected against partially embedded groups 
of younger stars (M43, ONC, OMC-2/3, L1641N). Image courtesy of Jon Christcnscn (christensenastroimages.com) ^ 
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Table 1. Catalogues and observations used in this study Table 2. CAHA 1.23m CCD observations 



Instrument 


Band / 




Channel 


XMM-Newton/EPIC 


0.1-10 keV 


SDSS 


u,g,r,i,z 


CFHT/MegaCam 


u,g,r 


2MASS 


J,H,Ks 


WISE 


3.3,4.6,12,22 ^m 


Spitzer/IRAC 


3.6,4.5,5.8,8.0 ^m 


Spitzer/MIPS 


24 fim 


Calar Alto 1.23m 


u,g,r 



Field 


RA (J2000) 


Dec (J2000) 


Trapezium 


nK.OK.in Q/ii 
0o:3t>:19.341 


— Uo:23:o0.oo 


Field 1 


05:35:24.651 


-05:55:06.69 


Field 2 


05:34:56.819 


-05:59:59.55 


Field 3 


05:35:25.044 


-05:59:15.44 


Field 4 


05:34:52.120 


-05:34:04.66 


Field 5 


05:33:59.980 


-05:35:40.17 


Field 6 


05:36:09.483 


-05:38:17.68 


Field 7 


05:37:23.116 


-05:56:10.97 



present a general discussion on the importance of the re- 
sult found in Sect. 4 and summarize the main results of the 
paper in Sect. 5. 



2. Data 

To characterize the foreground population to the Orion 
A molecular cloud we will make use of existing sur- 
veys together with raw data from Canada-France-Hawaii 
Telescope (CFHT), Calar Alto Observatory (CAHA 
1.23m), and the Spitzer satellite, that were processed and 
analyzed for the purpose of this investigation. 



2.1. Catalogues 

We retrieved the astrometry and photometry for all sources 
within a box of 5°xl5° centered around RA=85.7 and 
Dec=-4°(J2000) in the Sloan Digital Sky Su rvey III, the 



Wi de-fie ld Infrared Survey Explorer (WISE Cutri & et 



al 



201 2J, the Thir d XMM-Newton serendipitous source 
"d Watson et al.|2009J) and the 2MASS catalogues ( |Skrutskie 
et al.||2006[ ). Table [lj gives an overview of the properties of 



these catalogues . 



2.2. CFHT/Megacam 

A mosaic of 2x2 pointings covering 2°x2° centered on 
the Orion Nebula Cluster (ONC) was observed with 



CFHT/Megacam ( |Boulade etal 20030 with the Sloan ugr 
filters on 2005 February 14 (P.I. Cuillandre). Figure [2] gives 
an overview of the area covered by these observations. The 
conditions were photometric, as described in the Skyprobe 
database ( |Cuillandre et al.]|2004 ). Seeing was variable, os- 
cillating between 1-2" as measured in the images. A to- 
tal of 5 exposures of 150 s (u-band), 40 s (g-band), and 
40 s (r-band) each were obtained at each of the 4 positions. 
The observations were made in dither mode, with a jitter 
width of a few arcminutes at each position. This allows 
filling the CCD-to-CCD and position gaps and correcting 
for deviant pixels and cosmic ray events. The images were 
processed using the recommen ded Elixir reduction package 
( Magnier fc Cuillandre||2004[ ). A perture photometry was 
then extracted using SExtractor (Bertin & Arnouts 1996 1 
and the photometric zero-points in the SDSS system were 
derived by cross-matching with the SDSS catalogue. The 
CFHT/Megacam observations complement the SDSS data 
in one critical aspect: they provide data for regions around 
bright stars and nebulae, in particular the Orion Nebula 
region that is missing in the SDSS data. 



2.3. Calar Alto /1. 23m CCD Camera 

Selected pointings of the ONC (see Table [2] and Fig. 
were observed on 2011 December 15 with the Calar Alto 
CCD camera mounted on the 1.23m telescope (hereafter 
CAHA123). The CCD camera is a 2kx2k optical imager 
with a 17' field-of-view. The Sloan filters available at Calar 
Alto vignet the field and reduce it to a circular 11' diameter 
field-of-view. These observations are meant to complement 
the CFHT and SDSS observations below their saturation 
limits (at ugr «12 mag), and in the vicinity of bright satu- 
rated stars. Short exposures of 0.1 and 5.0 s were obtained 
in the Sloan gr filters, and of 0.1 and 10 s in the Sloan u fil- 
ter. The telescope was slightly defocused to avoid saturation 
of the brightest sta rs. Three standard fields (SA 97, SA 92 
and BD+21D0607, |Smith et al.|2002[ ) were observed during 
the course of the night to derive accurate zero-points. Each 
pointing was observed with a small dithering of a couple of 
arcminutes in order to correct for deviant pixels and cosmic 
ray events. The images were pre-processed (bias subtraction 
and flat-field correction) us ing standard pro cedure with the 
Eclipse reduction package (Devillard 1997). The astromet- 
ric registration and stacking were then per formed using the 
AstrOmatic software suite (Bertin 2010a). Aperture pho- 
tometry was finally extracted using SExtractor and the 
photometric zero-points in the SDSS system were derived 
by cross-matching with the SDSS and Megacam catalogues. 
The night was clear but not photometric. We observe a dis- 
persion in the zero-point measurements through the night 
of 0.06 mag in u and g, and 0.16 mag in r, which we add 
quadratically to the photometric measurement uncertain- 
ties. 



2.4. Spitzer IRAC 

The ONC has been extensively observed with IRAC on- 
board the Spitzer observatory in the course of programs 
30641, 43 and 50070. We retrieved the corresponding IRAC 
BCD images and associated ancillary products from the 
public archive, and processed them following stan dard pro- 
cedures with the r ecommended MOPEX software ( |Makovoz 



& Marleau |2005 ). The observations were all marie using 
the High Dynamics Range mode, providing short (0.6 s) 
and long (12 s) exposure. We processed the two sets inde- 
pendently so as to cover the largest dynamic range. The 
procedure within MOPEX includes overlap correction, re- 
sampling, interpolation (to have an output pixel scale of 
0'.'6) and outlier rejection. The individual fra mes were then 
median combined using Swarp (|Bertin|2010b|) using the rms 
maps produced by MOPEX as weight maps. Aperture pho- 
tometry of all the sources brighter than the 3-cr noise of the 
local background was extracted using SExtractor. We ver- 
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Fig. 3. g vs. g — r color-magnitude diagram for the entire survey in regions of increasing total line-of-sight extinction. 



ify that the corresponding ph otometry is in good agreement 
with IRAC photometry from Fang et al. ( 2009 ) within the 
uncertainties. 



2.5. Spitzer MIPS1 

The ONC was observed with the Spitzer Space Telescope 
and its MIPS instruments in the course of programs 202, 
30641, 30765, 3315, 40503, 47, 50070, and 58. We re- 
trieved from the public archive all the corresponding MIPS1 
(24 /zm) BCD images, and processed them with the rec- 
ommended MOPEX software. The procedure includes self- 
calibration (flat-fielding), overlap correction, outlier re- 
jection, and weighted coaddition into the final mosaic. 
Aperture photometry of all the sources brighter than the 
3-cr noise of the local background was extracted using 
SExtractor. We also verify that the corresponding photom- 
etry is in very goo d agreement with MIPS photometry from 
Fang et al. (2009) within the uncertainties. 



3. Results 

We are interested in the foreground population to the Orion 
A cloud, in particular the foreground populations towards 
the ONC. To separate it from the background we will use 
the optical properties of dust grains in the Orion A cloud 



to block the optical light to the cloud background. This 
is a very effective way of isolating the stellar population 
between Earth and the Orion A cloud, in particular if we 
use blue optical bands where dust extinction is most ef- 
fective. To select the final sample of foreground stars we 
take two filtering steps informed by Color-Magnitude and 
Color-Color diagrams in the optical and infrared. In par- 
ticular, we start by 1) using blue optical magnitudes and 
colors to define a reliable subsample of sources in front of 
the cloud, then 2) using a near-infrared color-color diagram 
to reject sources affected by extinction (these are sources 
that are either young stars inside the cloud or background 
sources that are bright enough to be detected in the optical 
survey) . 

The first filtering step is displayed in Figure [3] where 
we present a g vs g — r color-magnitude (CMD) diagram 
combining the SDSS, CFHT/Megacam and CAHA123 pho- 
tometry. We chose g — r over the more extinction sensitive 
u — g as the u— band observations are significantly less deep, 
and extremely sensitive to excesses related to accretion and 
activity. 

The first diagram (on the left) represents the g vs g — r 
color-magnitude diagram for all the sources in our com- 
bined database. The three CMD-diagrams to the right are a 
subset of the first containing only sources projected against 
increasing contours of dust extinction of the Orion A cloud 
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Fig. 4. J— H versus H 
giant s equences from 
(20071. The sizes oft 



-K Color-Color vers us J brightness diagram. The solid grey lines represent t he main sequence and 



Bessell fc Brett ( 1988 ) . The dashed grey line represents the main sequence from Kraus & Hillenbrand 
le symbols are proportional to J-band brightness. Sources taken as foreground candidates are marked 



in blue, while rejected sources, mostly extincted sources, are marked in red. 



(about 4, 5, and 6 magnitudes of visual extinction). These 
colum n density threshold s were estimated from the 13 CO 
map of Bally et al. ( 1987[), cross-cali brated with the extinc- 
tion map of Lombardi et al. (|2011[). W hile using directly 



the extinction map of 
lar results, we preferre* 



(2011) gives simi- 



^ jLombardi et al. 

d to avoid dealing with any possible 
systematics affecting this map caused by a potential sub- 
stantial population of foreground sources. As we impose the 
condition of keeping only sources that are seen against in- 
creasing levels of dust extinction two things occur: 1) the 
number of stars decreases, naturally, because the solid an- 
gle on the sky decreases, and 2) a well defined sequence 
appears. This sequence is not what is expected from the 
general Galactic population between Earth and the Orion 
A cloud at 400 pc, as confirmed with the Besacpn stellar 
population model A Robin et al. 



retain the subsample of sources 



2003) 
that 



From this step we 
is seen in projection 



against column densities of greater than ~ 5 visual mag- 
nitudes of extinction (third panel in Figure [3]) , or a total 
of 2169 sources from more than 1.25xl0 5 sources in the 
combined SDSS-MEGACAM-CAHA123 catalog. Most of 
the discarded sources have colors consistent with unrcd- 
dened and slightly reddened unrelated field stars towards 
the background of the Orion A cloud. Among the sources 
that pass the first filter there could be some with g-band 
excess emission, but these should have a negligible effect in 
the selection process, in particular because the next filter- 
ing step is done at the near-infrared. 

The second filtering step consists of discarding extincted 
sources. We want to remove from the sample any source 
that might be associated with the cloud (young stellar ob- 
jects still embedded in the cloud, for example), as well as 
background sources that are bright enough to be detected 
at ~0.487/im (g-band) through Ay ~ 5 mag of cloud mate- 
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Foreground sources Reddened sources 




Fig. 5. Left panel: spatial density of the foreground sources (blue sample). The unshaded area in the Figure represents 
the region of the cloud where Ay > 5 mag, on which the selection was performed. The blue contours (with increments of 
10% from the maximum) represent the surface density of foreground sources (constructed with a gaussian kernel with a 
width of 20'). Right panel: same as in the left panel but for the reddened sources. The distribution of foreground sources 
shows a well defined peak coinciding with the poorly studied NGC 1980. The reddened sources, on the other hand, peak 
around 1) the Trapezium cluster and are mostly confined to the nebula and 2) the L1641N star forming region, with a 
peak towards a hitherto unrecognized group of YSOs (see text). The reddened and foreground populations are spatially 
uncorrelated but there is significant overlap between the two, in particular with the sources inside the Orion Nebula. 



rial. We perform this filtering using a J— H vs. H— K Color 
Color diagram where extincted sources are easily identi 
fied along the reddening band, away f rom their unreddened 
main sequence (and giant) colors (e.g.jLada k Adams|1 992 



Alves et al.||1998l iLombardi k Alves|200ip . We present in 
Figure |4j the J— H vs. H— K Color-Color diagram for the 
2169 sources that passed the first filtering step. The size 
of the symbols in this Figure are proportional to J-band 
brightness. The selection criteria used to identify the likely 
foreground population was: 



„ 0.96 -(J-H) 
H < ■ r -^ z mag 



1.05 



J < 15 mag 



(1) 



(2) 



J-H < 0.74 U H - K > -0.2 U H - K < 0.43 mag (3) 

Sources that are consistent with having no extinction 
within the photometric errors, are marked in blue, while 



rejected sources are marked in red. Condition (1), the main 
filter, is taken as the border between extincted and non 
extincted sources, and it was selected to be roughly par- 
allel to the main-sequence early M-star branch (to about 
the color of a M4-M5 star). Condition (2) and (3) further 
reject sources that are faint or have dubious NIR colors 
(either bluer than physically possible, or redder than main- 
sequence stars, suggestive of a NIR excess). Condition (2) 
in particular makes the selection more robust against pho- 
tometric errors (the typical photometric error imposed by 
condition (2) is J err ~ 0.03±0.01, R err ~ 0.03±0.01, and 
K err ~ 0.03 ± 0.01 mag, which translates into a maximum 
Ay error of ~ 1 mag), and should reach a sensitivity limit 
capable of including M4 main-sequence stars at the distance 
of the cloud (J ~ 15 mag). More than two thirds of the 2169 
sources are rejected (red sources) and a total of 624 sources 
have colors consistent with foreground stars suffering no or 
negligible amounts of extinction)^] 



1 A table of candidate foreground sources is provided in elec- 
tronic format. 
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X-ray 

(any XMM- Newton detection) 




Fig. 6. Spatial density of all the X-ray sources in the Third 
XMM-Newton serendipitous source catalog for the same 
region as in the previous Figure. The gaussian kernel used 
and contour separation are the same. NGC 1980 is detected 
as a distinct enhancement in the surface density of X-ray 
sources, together with the Trapezium cluster, the L1641N 
population, and the hitherto unrecognized group that we 
name L1641W. 



3.1. Spatial distribution of the foreground and reddened 
sources 

The addition of the third dimension ( J-band brightness) to 
Figure [4] makes obvious the presence of a well populated 
main-sequence branch, from early types (B-stars) to late 
types (M-stars). The clear presence of so many early- type 
stars, as well as a well defined main sequence, suggests that 
the foreground population is not dominated by the random 
Galactic field from Earth to the ONC region, but could in- 
stead be a well defined stellar population. To investigate 
this idea we present in the left panel of Figure [5] the spatial 
density of foreground sources (blue sample) . The unshaded 
area in the Figure represents the region of the cloud where 
Ay > 5 mag, on which the selection was performed. We 
constructed a surface density map of foreground sources 
using a gaussian kernel with a width of 20' represented in 
the Figure as blue contours (with increments of 10% from 
the maximum). It is clear from Figure p\ that the distri- 
bution of foreground sources is not uniform, as expected 
for the Galactic field between Earth and the ONC region 
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Fig. 7. a) K-band luminosity functions for the foreground 
(blue) and the extincted (red) sample (see Figure [4]). Albeit 
being affected by extinction, the red sample is surpris- 
ingly brighter than the foreground sample (blue), suggest- 
ing that the two populations are intrinsically different, b) 
Empirical cumulative distribution functions for both sam- 
ples, together with upper and lower simultaneous 95% con- 
fidence curves, confirming that the two populations have 
statistically different luminosity functions. 



in the Orion A cloud, but is instead strongly peaked and 
fairly symmetric. The peak of the distribution coincides 
spatially with the poorly studied iota Ori cluster, or NGC 
1980, suggesting that the foreground population is domi- 
nated by NGC 1980 cluster members. The elongated shape 
of the peak is not meaningful as it is caused by the rela- 
tive narrow Ay > 5 mag region on which the density was 
calculated. 
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The right panel of Figure [5] shows the spatial distribu- 
tion of reddened sources. This distribution is dominated by 
two peaks, a relatively well defined one in the Trapezium 
cluster region and another, more diffuse, coinciding with 
the overall gas distribution around the L1641N star form- 
ing region. This is not surprising as the reddened sources 
are expected to be dominated by the embedded young stel- 
lar objects in the Orion A cloud, and these are know to 
cluster around these two regions. More striking, instead, is 
that the foreground sources and the reddened sources ap- 
pear spatially anti-correlated: the maximum of the distribu- 
tion of the foreground sources coincides with the minimum 
of the distribution of reddened sources. This is evidence 
that the foreground population is not the emerging young 
stellar population from the ONC but is instead an entirely 
different population. Because it contains a fully sampled 
and unreddened main-sequence from B- to M-stars (see 
Figure |4|, this population is most likely the stellar pop- 
ulation of the NGC 1980 cluster, seen in projection against 
the Orion A cloud. But because this population overlaps 
significantly with the ONC (the distribution of foreground 
sources appears symmetric to about 7 pc from its center), 
it implies that the ONC is not comprised by a single cluster 
with the Trapezium as its core, but has instead three stel- 
lar populations (a) the youngest population, including the 
Trapezium and ongoing star formation associated with the 
dense gas in the nebula, b) part of the NGC 1980 cluster in 
its foreground, and c) the unrelated, Galactic, foreground 
and background population. 

In Figure [6 we present, for the same region shown in 
Figure [5] the distribution of all the X-ray sources in the 
Third XMM-Newton serendipitous source catalog. We want 
to investigate if the NGC 1980 X-ray source counts ap- 
pear as a distinct peak in the surface density map (con- 
structed using a gaussian kernel with a width of 20' as in 
the previous Figure). As can be seen in this Figure, both 
the foreground and the extincted populations are detected, 
and a distinct enhancement in the surface density of X- 
ray sources is seen towards the center of NGC 1980, given 
strength to the idea that the foreground population is not 
the emerging young stellar population from the ONC but is 
instead an entirely different population. The highest peak 
in this surface density map is centered on the Trapezium 
cluster. This was not the case for the extincted population 
as seen in Figure [5j which peaks slightly to the South of 
the Trapezium, although this mismatch could be due sim- 
ply to the fact that our MEGACAM g-band observations 
are affected by the bright nebula and dust extinction in the 
Trapezium region, hence substantially less sensitive to the 
embedded population of the Trapezium cluster. 

3.1.1. A hitherto unrecognized group of YSOs? 

The presence of an enhancement of the reddened sources 
in Figure [H] towards RA: 5 h 35 m , Dec: -6°18 m , immedi- 
ately South of NGC 1980 and towards the West of what is 
norm ally taken as the L1641N cluster (e.g. Allen & Davis 
2008 1, is tantalizing. Could this relatively small enhance 



ment be another hitherto unrecognized group of YSOs? The 
enhancement is clearly detected in X-rays (Figure [6]), and 
is tentatively detected in the optical (Figure pi), providing 
support in favor of this possibility. We name this poten- 
tially new group of about 50 stars (counted on the reddened 
sample) as L1641W. The group is not associated with any 



Table 3. Position of clusters in Figure [(jj including the 
newly identified L1641W 



Field 


RA (J2000) 


Dec (J2000) 


L1641W 


05:34:51.0 


-06:17:40 


NGC 1980 


05:35:11.0 


-05:58:00 


Trapezium 


05:35:16.5 


-05:23:14 


L1641N 


05:35:55.7 


-06:23:55 



obvious nebula nor does it include any obvious bright star. 
Because it appears less extincted than the L1641N popu- 
lation, a nd it is not obviously detected in the Spitzer sur- 
vey (e.g. I Allen fc Davis||2008 ), suggests that it is probably 
more evolved than the L1641N population. We speculate 
that this new group is either a foreground young group 
ramming into the Orion A cloud, or a slightly older sibling 
of NGC 1641N, leaving the cloud. 



3.2. Luminosity function of foreground and reddened sources 

In Figure vf\ a) we present the K-band luminosity functions 
for both the foreground (blue) and reddened (red) samples 
(see FigureB. To enable a direct comparison, the reddened 
sample was also constrained with condition (2), namely, J < 
15 mag. Surprisingly, the extincted sample (red) is brighter 
than the foreground sample (blue), even if no derredening 
procedure was applied to the red sample. We confirm that 
the differences between the two luminosity functions are 
significant, to a 95% confidence level, by analyzing their 
empirical cumulative distribution functions (see Figure [7] 
b)). Note that had we de-reddened the extincted sample, 
the difference between the luminosity functions would have 
been even higher. This suggests, like in the previous sec- 
tion, that the foreground and reddened population are in- 
trinsically different. A likely explanation for the difference 
in the luminosity functions is that the reddened sample is 
dominated by very young stellar objects still embedded the 
cloud, which are intrinsically brighter than normal stars be- 
cause of both stellar evolution and the presence of K-band 
excess emission. 



3.3. Velocity dispersion profile 



Tobin et al. (2009) presented a large kinematic study of 
the OWC, covering about 2° of Declination centered on the 
Orion A cloud, from NGC 1977 down to L1641N. This sur- 

and 



(2008) 



vey builds on previous work by Furesz et al. 
constitutes the largest and highest precision kinematic sur- 
vey of this region to date, offering a unique possibility to 
characterize kinematically the ONC foreground population 
identified in this paper. In Figure [8] we present the North- 
South velocity dispersion profile of the ONC region, taken 
from Table 13 of [Tobin et al] ( |2009[ ). The filled circles repre- 
sent the sources m this paper with reliable radial velocities. 
The NGC 1977, the OMC2/3, the Trapezium cluster, and 
NGC 1980 are indicated, as well as the extent of the Orion 
Nebula (light open circle). The thick red line represents the 
North-South velocity dispersion profile measured in bins of 
Declination (indicated by the thin horizontal lines). 

It is striking that the velocity dispersion profile has a 
minimum at the location of NGC 1980. This is perhaps the 
strongest indication we have that the stellar population of 
NGC 1980 is a distinct population from the reddened pop- 
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Fig. 8. The North-South velocity dispersion profile of the 
ONC region . The filled circles represent the sources in |Tobin| 
et al.| ( |2009| ) with reliable radial velocities. The NGC 1977, 
the Trapezium cluster, and NGC 1980 are indicated, as well 
as the extent of the Orion Nebula (light open circle). The 
thick red line represents the North-South velocity disper- 
sion profile measured in bins of Declination (indicated by 
the thin horizontal lines). There is an increase from the 
Trapezium to the edge of the nebula to the South, followed 
by a clear minimum around NGC 1980, strongly suggesting 
that NGC 1980 is a well defined and different population 
from the st ellar population in side the nebula. Data from 
Table 13 oflTobin et all (12009 1. 



ulation inside Orion A. The measurement of velocity dis- 
persion in the bin that mostly includes NGC 1980 (a =2.1 
km/s) was not optimized to isolate the most probable mem- 
bers of this cluster, and should then be seen as an upper 
limit to the true velocity dispersion in this cluster. Still, this 
value is close to the velocity dispersion of the Trapezium 
cluster as measured from the proper motion of stars within 



one half degree of the center of the Trapeziu m, namely, 
1.34±0.18 km/s for a sample of brighter stars (van Altena 
et al.|198"8| and 1.99±0.08 k m/s for a larger sample includ- 
ing relatively fainter stars (Jones & Walker [1988 ). Both 
velocity dispersions were corrected to the more recently es- 
timate of the distance to the Orion A cloud (400 pc) . 

Given the striking differences in the velocity dispersion 
profile we then calculated the mean radial velocity per bin 
from the subsample of single sources (not directly available 
in the Tobin et al. 2009 paper) and found that although 
showing variations from bin to bin, these variations are of 
the order of the measured dispersions. In particular, the 
mean velocity for the bin including the Trapezium (—5.3° < 
5 < -5.4°) and NGC 1980 (-5.8° < 6 < -6.0°) is 25.7±3.0 
km/s and 24.3 ± 2.7 km/s respectively. Within the errors, 
estimated as the median absolute deviation in each bin, the 
NGC 1980 cluster has virtually the same radial velocity as 
the ONC. We note, however, that we are taking the bins 
as simple slices at constant declination, without trying to 
optimize their boundaries to better separate the different 
populations. 

Because of the importance of measuring the velocity dif- 
ferences between the Trapezium and NGC 1980, especially 
for a discussion on the origin of NGC 1980, we made an al- 
ternative source selection and created two new subsamples, 
that are in principle more pure, but have about three times 
less sources. For the NGC 1980 subsample we matched the 
Tobin et al. 2009 catalog with the foreground population 
identified in this paper. For the Trapezium we m atched the 
Tobin et al . 2009 catalog with the COUP sample ( Feigelson 
et al.||2"002" ), that is dominated by Trapezium sources, and 
removed sources that matched the foreground population. 
Because this Trapezium subsample should be of "high con- 
fidence" , we used the radial velocity limits found in this 
subsample (6.2 km/s and 36.6 km/s) to exclude 5 extreme 
outliers in the NGC 1980 sample (with velocities of ~ —40 
and ~ 90 km/s). In these subsamples, the mean velocity for 
the Trapezium and NGC 1980 clusters is 25.4 ± 3.0 km/s 
and 24.4 ± 1.5 km/s respectively, or essentially the same 
values as derived above, with the important difference that 
the dispersion of velocities in NGC 1980 is now reduced 
by about a factor of two, once again suggesting that this 
cluster is a distinct population from the reddened popu- 
lation inside Orion, as argued above. Still even with the 
decreased velocity dispersion the measured velocity differ- 
ence of 1 km/s is not statistically significant. 

3.4. On the age and population size of NGC 1980 

In order to estimate an age to the NGC 1980 cluster 
we compare the evolutionary status of class II sources in 
various clusters through the analysis of the median spec- 
tral energy distribution (SED) of late-ty pe (spectral type 
later t han K0) members. We follow the Hartmann et al. 
(2005) definition of Class II, namely objects with 0.2 < 
[4.5] < 0.7 mag and 0.6 < [5.8] - [8.0] < 1.1 mag. 
To compute the median SED for the different clusters 
we retrieved the optical, near-infrared (2MASS) and mid- 
infrared (Spitzer and WISE) photome try for samples of con 
firmed members of Taurus (1-3 Myr, 
IC 348 (1-3 Myr, |Lada etaT " 

„i ll 'om ' AI l /"l„4. 



Winston et al. 2010| Hutermut 



2006 



(1 



Myr, 
Ori (5- 

7 My r, |Dolan fc Mathieu||2001| |Barrado y Navascues et al 



i et al. 



2004), and rj-Cha (5-10 Myr, 



Luhman et al. 2010) 



NGC1333 
2008] 



A- 



Megeath et al.||2005[ ). To 
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compute the median SED for the Trapezium cluster we de- 
fined first a "high confidence" Trapezium member catalog, 
as we did in the previous Section, by cross-mat ching the 
X-ray COUP sample from (Prisinzano et al. 2008) with the 
foreground (NGC 1980) sample, and excluding all matches 
as unrelated foregrounds. The individual SEDs within each 
cluster were normalized to the J-band flux, and the median 
cluster SED of each cluster was computed. Figure [9] shows 
the result. One can see from this Figure [9] that the optical 
part of the SED varies from cluster to cluster, mostly due to 
dust extinction. More striking, the mid-infrared (> 3 ^m) 
excesses, related to the presence of a disk, decrease system- 
atically with age. 

The median SED of NGC 1980 seems to fit between the 
median SED of Taurus (1-3 Myr) and A-Ori (5-7 Myr), 
suggesting an age in between that of these regions. But 
another constraint is given by the massive stars in the cen- 
ter of the cluster. Of the five brightest stars at the peak 
of the spatial distribution in Figure [5] only the brightest, 
iota Ori (09 III, V=2.77 mag), seems to have evolved from 
the main sequence. This implies an age of about 4-5 Myr 
for t his star, assum ing it started its life as a 25 M star 



(e.g. Massey 2003). This age fits well within the inferred 



age from the median SED and is also in agreement with 



the estimate of Warren & Hesser ( 1978 ) for the age of Ori 
OB lc subgroup (of about 4 Myr) . 

To estimate the size of the cluster population we concen- 
trate on the distribution of foreground sources from the cen- 
ter of the cluster to the South in order to avoid incomplete- 
ness issues caused by the bright Orion Nebula. We counted 
the number of sources falling on a 20° "pie slice" inside the 
Ay > 5 mag region, centered on the cluster and with a ra- 
dius of 7 pc. This radius corresponds approximately to the 
extent of the 10% contour in Figure[5] chosen to account for 
contamination from the Galactic field between Earth and 
Orion (est imate d to be 6-9% of the foreground population 
in Section 3.5.1). Note that this radius is not the half-mass 



radius but simply the radius to which we can trace the en- 
hancement of sources over the unrelated foreground field. 
We repeated this measurement several times to account for 
uncertainties in the location of the cluster center and ob- 
tained an average number of 110 sources in the 20° "pie 
slice" . Assuming spherical symmetry for the distribution of 
sources in NGC 1980, we expect then a total of about 2000 
sources in NGC 1980, or a total cluster mass of about 1000 
M Q (assuming an average mass per star of 0.5 M Q ). 

Assuming NGC 1980 has a normal Initial Mass Function 
(IMF), we can make a consistency test on the likelihood of 
the number of sources in this cluster being of the order of 
2000. For this we constructed 200000 synthetic clusters of 
1000 M each by randomly sampling the Kro upa and the 
Chabrier IMFs ( |Kroupa|[200l| |Chabrierl 



2003 ) and tracked 



the mass of the most massive star in each synthetic cluster. 
The mean mass of the most massive star was 54±26 M 
(Kroupa) and 22±11 M (Chabrier). Assuming there were 
no supernovae in NGC 1980 yet, and that iota Ori is the 
most massive star in the cluster then, to first approxima- 
tion, a population of about 2000 sources seems a plausible 
estimate of the size of NGC 1980 population. 

At the moment we cannot derive a reliable cluster radial 
profile, nor a half-mass radius, or even be certain about the 
position of the center of the cluster, as our optical observa- 
tions are incomplete in the vicinity of the early type stars of 
NGC 1980, and we only have a "pie slice" view on the radial 
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Fig. 9. Median SED of class II sources in clusters of 
various evolutionary stages: NGC1333 (magenta, 1 Myr), 
Trapezium COUP sources (black, 1 Myr), Taurus (green, 
1-3 Myr), NGC1980 members (red), A Ori (blue 5-7 Myr) 
and 77— Cha (orange, 5-10 Myr). 



extent of the cluster. This should be improved in follow-up 
work, in particular in combination of dedicated NIR ob- 
servations which are less sensitive to the large brightness 
contrasts between early and late type stars in this cluster. 



3.4.1. On the origin of NGC 1980 and impact on the Orion A 
cloud 



We found in Section 



3.3| that the radial velocity of NGC 
1980 is indistinguishable, or has a difference of the order of 
a few km/s at best from the radial velocity of the embed- 
ded Trapezium population. This surprising result implies 
that the radial velocity of NGC 1980 is essentially the same 
as the velocity of the gas in the Orion A cloud, since the 
ONC population h as the same radial velocity as the cloud 
( |Tobin et al.1|2009[ ). This strongly suggests that NGC 1980 
is somehow connected to the Orion A cloud, or better, that 
the cloud that formed NGC 1980 was physically related to 
the current Orion A cloud. One would not expect the dis- 
tance to NGC 1980 to be substantially different than the 
current distance estimate to the ONC, and a fitting of the 
ZAMS on the optical data presented in this paper is indeed 
consistent with a distance of 400 pc. 

Despite its relatively older age, lack of obvious H n re- 
gion, and lack of measurable dust extinction, NGC 1980 
moves away from Earth at the same velocity as the large 
Orion A cloud on which it is seen in projection. Because of 
their likely proximity, one wonders on the effects of the ion- 
izing stars from NGC 1980 on the Orion A cloud, or what 
the cloud was about 4-5 Myr ago, in particular on a possible 
acceleration and compression of the cloud by the UV radi- 
ation from these stars. How important was/is this process 
in this region? Could the formation of the O NC have bee n 
triggered by its older sibling, as suggested in Bally (2008)? 



At first look our results would argue that the impact would 
have been minimal, NGC 1980 has essentially the same ra- 



dial v el ocity as the Orion A cloud, but the work of |Oort 
(|1954[); |Kahn| (|1954[); |Oort & Spitzer| (|1955[) suggests that 
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Fig. 10. Matches bet ween the foreground contamination and well known catalogs of the ONC. a) In green, the optical 
spectroscopic survey of|Tobin et al.| J2009 ), with the matches with the foreground population marked in blue, b) In yellow, 



the well studied sample of Hillenbrand (19971), with the matches with the foreg round population marked in blue, c) In 

l theCi 



red, the sample of X-ray sources from the COUP project (Feigelson et al. [2002 1 , with the matches with the foreground 
population marked in blue. As in Figure [5] the unshaded area represents the region of the cloud where Ay > 5 mag, on 
which the selection of foreground sources was performed. 



final speeds between the ionizing star and the cloud would 
be of the order of a few km/s, which cannot be ruled out by 
the current accuracy of the data. While our results do not 
give final evidence in support of the tantalizing suggestion 
that the formation of the ONC could have been triggered 
by NGC 1980, they are not inconsistent with it either. A 
new dedicated radial velocity survey of the region, together 
with a sensitive proper motion survey, are needed to un- 
derstand the interplay between these two massive clusters. 
This configuration (an embedded cluster in the vicinity of 
a ~ 5 Myr cluster) is unlikely to be unique in massive star 
forming clouds, but it will best addressed in the nearest 
example. 



3.5. Contamination of ONC catalogs 

We have showed above that there is a rich and distinct 
foreground population of stars, likely associated with the 
young (~ 5 Myr) poorly studied but massive NGC 1980 
cluster, that is not directly associated with the ongoing 
star formation in the ONC. This finding raises concerns 
on the contamination of currently available observables for 
this important region, and future studies should take this 
foreground population into account. But how large is this 
contamination? There are two we ll known ONC cat alogs 
used in the lit erature, namely the [Hillenbrand| (|1997l) and 



the catalogs of|Da Rio et al.| ( |2009| ); pa Rio et al.|P0T0| );|Da 
Rio et al. (2012) covering a roughly square area of about 
0.5" x 0.5 (~ 3.5 x 3.5 pc) cen tered on the Trapezium 
cluster. The Da Rio et al. (20121 supersedes all previous 
catalogs, but it is the most recent hence less used in the 



alog has been used extensively in the literature and has 
spawned a large number of the star formation studies on 
the star formation properties of the ONC region. We esti- 
ma te here the likely fo reground contamination fraction for 
the Hillenbrand (1997) catalog as it is the most used one, 
but also because it is likely to be the least contaminated 
since the Da Rio catalogs cover a slightly larger area of the 
sky towards NGC 1980. 

To estimate the probable contamination fraction of 



Hillenbrand (1997) we matched the foreground population 
with this catalog for stars falling within the Ay > 5 mag 
region where the foreground was selected (see Figure pj 
and where I-band < 16 mag. The last constraint accounts 
for the fact that the Hillenbrand ( 1997 ) sample is not uni- 
formly deep (it reaches about 2 magnitudes deeper around 
the Trapezium cluster) , and that the selection of foreground 
stars, made at g-band, seems complete to about I-band ~ 
16 mag (af ter transformation o f the SDSS photometry into 
Johnson's fllvezic et all 120071)). We find that 11% of the 
' feTHill - 



sources m 



brand 



(1997) catalog have a match in 
the foreground sample (8% if we remove the constraint on 
the I-band brightness). If one sees the Trapezium cluster as 
a component of the ONC, and not as the only component, 
and remove it from consideration, then the fraction of fore- 
ground contaminants in the ONC rises to 32%. For this esti- 
mate the area on the sky covered by the Trapezium cluster 
is taken as an ellipse with a= 7.5' and b= 3.8', with a posi- 
tion angle of —10°, s imilar to the definition in Hillenbrand 
|fc Hartmann] ( |1998[ ). One can also make an estimate of 



the possible contamination to the entire Hillenbr and] (|1997[) 
catalog by applying equation (jTj) to the entire Hillenbrand 



community. On the other hand, the Hillenbrand (1997) cat- (1997) catalog, which gives then a contamination fraction 
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of 20%, or 63% when the Trapezium is removed from con- 
sideration. Note that all these estimations assume that the 
fraction of ONC stars without measurable extinction is neg- 
ligible, which is likely given the distribution of foreground 
stars in Figure [5j but will need to be investigated further 
in future work. 

Even without removing the Trapezium cluster from con- 
sideration, contamination fractions of the order 10-20% are 
significant and will necessarily lead to systematic errors in 
the basic derived physical quantities for this star forma- 
tion benchmark. Still, these are necessarily lower limits to 
the true contamination fraction of the ONC sample for at 
least two reasons: 1) our g-band MEG AC AM survey is not 
as sensitive in regions of high nebula brightness, especially 
around the Trapezium, and 2) we are not sensitive, by de- 
sign, to background sources. While it is normally argued 
that the high background extinction behind the Trapezium 
blocks most background stars, this is only valid for the inner 
regions of the Trapezium cluster (~ 25' 2 ), but not valid for 
the entire ~ 700' 2 Orion Nebula, (e.g. Johnstone & Bally 
1999| jBuckle et al.|2012|> So back ground contamination is 
variable across the ON C and expected for any optical or 
infrared survey of this region. In regards to the contamina- 
tion of the ONC region by NGC 1980, it is a function of the 
position in the Nebula, having a minimum at the center of 
the Trapezium where the Trapezium cluster stellar density 
is highest, increasing gradually towards the South as one 
approaches the core of NGC 1980 (see Figure [§. 

3.5.1. Unrelated Galactic field foreground population 

Most of the identified foreground population (624 sources) 
is likely to belong to NGC 1980, as seen from the symmet- 
ric and peaked spatial distribution in Figure [5| but some 
fraction of these must be made of the Galactic field pop- 
ulation between Earth and the Orion A cloud. A first and 
simple estimate of the size of this population can be made 
by correlating the f oreground population with the sample of 



Tobin et al. (120091) for which good radial velocity measure- 
ments exist (see Figure 10). The distribution of radial ve- 



locities of the 188 sources m common reveals a gaussian like 
distribution centered at ~ 26.1 km/s, with a gap roughly 
between —20-0 km/s and 40-60 km/s without any source, 
and 3 sources below —20 km/s and 9 above 60 km/s (i.e., 
12 potential outliers). Not surprisingly, if we sigma-clip the 
entire distribution at 3 er we find 11 outliers. If we remove 
from the distribution the 12 potential outliers described 
above and then sigma-clip the rest of the distribution we 
find 5 more outliers, but this time at the wings of the dis- 
tribution. This suggests that about 6-9% of the foreground 
sources identified in this work are likely field sources un- 
related to NGC 1980. This estimate of the Galactic field 
foreground contamination for Orion A is in rough agree- 
ment with what would be expected from the Besacon stel- 
lar population model (Robin et al. 2003) for the depth of 
our MEGACAM survey. 

4. Conclusions 

We have made a link between the foreground population 
towards the well known star formation benchmark ONC re- 
gion and the stellar population of the poorly studied NGC 
1980 cluster (or iota Ori cluster). Not only did we detect the 
presence of a well populated main-sequence (from B-stars 



to M-stars), the foreground sources have: 1) a well defined 
spatial distribution peaking near iota Ori, 2) a fainter lu- 
minosity function when compared to the extincted young 
population embedded inside the cloud, and 3) a low velocity 
dispersion, typical of that of other young clusters. 

Unlike the ONC, NGC 1980 is a relatively older clus- 
ter (4 — 5 Myr), lacks an obvious Hn region, and is com- 
paratively free of dust extinction. Surprisingly, the radial 
velocity of NGC 1980 is currently indistinguishable from 
the radial velocity of the ONC embedded population or the 
radial velocity of Orion A cloud, suggesting that both clus- 
ters are genetically related, and at about the same distance 
from Earth. 

A general concern that this study raises is the risk of 
population mixing in star formation studies. It is unlikely 
that the ONC is atypical in this respect and a dedicated 
multi-wavelength study to disentangle the different popu- 
lations, together with a sensitive proper motion survey of 
the region, is urgently needed. The ONC is still the closer 
massive star formation region to Earth, and albeit more 
complicated than first assumed, it is still the one offering 
the best detailed view on the formation of massive stars 
and clusters. 

The main conclusions of this work are: 



— We make use of the optical effects of dust extinction to 
block the background stellar population to the Orion 
A cloud, and find that there is a rich foreground stel- 
lar population in front of the cloud, in particular the 
ONC. This population contains a well populated main- 
sequence, from B-stars to M-stars. 

— The spatial distribution of the foreground population 
is not random but clusters strongly around NGC 1980 
(iota Ori), has a fainter luminosity function, and dif- 
ferent velocity dispersion from the reddened population 
inside the Orion A cloud. This foreground population 
is, in all likelihood, the extended stellar content of the 
poorly studied NGC 1980 cluster. 

— We estimate the number of members of NGC 1980 to 
be of the order of 2000, which makes it one of the most 
massive clusters in the entire Orion complex, and es- 
timate its age to be ~ 4—5 Myr by making a com- 
parative study of median spectral energy distributions 
among known young populations and constraints from 
the age of post main sequence star iota Ori. 

— This newly found population overlaps significantly with 
what is currently assumed to be the ONC and the 
L1641N populations, and can make up for more than 
10-20% of what is currently taken as the ONC popula- 
tion (30-60% if the Trapezium cluster is removed from 
consideration). 

— Our results suggest that what is normally taken in the 
literature as the ONC should be seen as a mix of sev- 
eral unrelated populations: 1) the youngest population, 
including the Trapezium cluster and ongoing star for- 
mation in the dense gas inside the nebula, 2) the young 
foreground population, dominated by the NGC 1980 
cluster, and 3) the poorly constrained population of 
foreground and background Galactic field stars. 

— We re-determine the mean radial velocity for the 
Trapezium and NGC 1980 clusters to be 25.4±3.0 km/s 
and 24.4 ± 1.5 km/s respectively, or indistinguishable 
within the errors, and similar to the radial velocity of 
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the Orion A cloud, suggestive of a genetical connection 
between the two. 

— We identify a hitherto unrecognized group of about 50 
YSOs West of L1641N (L1654W) that we speculate is 
cither a foreground group ramming into the Orion A 
cloud, or a slightly older sibling of NGC 1641N, leaving 
the cloud. 

— This work supports a scenario where the ONC and 
L1641N are not directly associated with NGC 1980, 
i.e., they are not the same population emerging from 
its parental cloud but are instead distinct overlapping 
populations. This calls for a revision of most of the ob- 
servables in the benchmark ONC region (e.g., ages, age 
spread, mass function, disk frequency, etc.). 
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